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Abstract 

We show that a class of Affleck-Dine baryogenesis directly relates the observed 
mass density of baryons, Ob, to that of dark matter, Q*dm- In this scenario, the 
ratio of baryon to dark matter mass density is solely determined by the low en- 
ergy parameters, except for an O(0. 1) effective CP-violating phase. We find that 
Qb/^dm = 0(0.1) with reasonable parameters, which lies surprisingly just in the 
range of observation. This scenario is totally free from the cosmological gravitino 
problem, and independent of the detailed history of the Universe as long as it sat- 
isfies quite weak constraints. 



On leave from University of Tokyo. 



1 Introduction 



The searches for both the origin of the observed baryon asymmetry and that of the 
dark matter in the present Universe have been attractive subjects for many physicists. 
The great success of the Big Bang nucleosynthesis (BBN) predicts an abundance of the 
baryon asymmetry in the present Universe as QbIi 2 ~ 0.02. Here, h is the present Hubble 
parameter in units of 100 km sec -1 Mpc -1 , and = Pb/Pc- (pb and p c are the energy 
density of the baryon and the critical energy density in the present Universe.) As for the 
dark matter, various observations of the dynamics of galaxies and their clusters, as well as 
theoretical analyses of the large-scale structure formation indicate that its cosmological 
abundance lies in the range, Qj)Mh 2 — 0.1-0.2. 

There exist a number of scenarios of "baryogenesis" , which can generate the required 
baryon asymmetry. Also for dark matter, several particles have been proposed as candi- 
dates, such as axion, neutrino and, among others, the lightest supersymmetric (SUSY) 
particle (LSP) (especially the lightest neutralino, \)- 

However, there remains one prominent problem. The matter density of baryons is 
roughly O(10-20)% of that of dark matter, i.e. Q b /^dm — 0.1-0.2. Why are they so 
close to each other? Most of the models of baryogenesis use baryon B- and/or lepton L- 
number-violating operators in high energy physics. On the other hand, the abundance of 
dark matter does not seem to be at all related with such high energy physics relevant to 
baryogenesis. In the standard SUSY dark matter scenario, for example, the mass density 
of the dark matter is solely determined by the weak scale properties of the neutralino LSP, 
such as the mass and the annihilation cross section. As long as dark matter is composed 
of the thermal relics of the LSPs, it is completely independent of other details in the 
history of the Universe. It must be very exciting if we can find a simple mechanism to 
directly relate the two independent physical parameters, the baryon asymmetry and the 
dark matter density. 

One interesting solution to this puzzle was proposed by using the Affleck-Dine (AD) 
baryogenesis [1[] in the context of the minimal supergravity (mSUGRA) scenario [0. [] 
In this scenario, non-topological solitons called Q-balls Q play the role of the common 
source of the baryon asymmetry and the LSP. In the AD baryogenesis, a complex scalar 

1 The possibility to explain both the baryon asymmetry and dark matter by a single mechanism was 
first considered in the context of gauge- mediated SUSY breaking scenarios M . 
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field 0, which is a linear combination of squark fields, obtains a large expectation value 
along one of the flat directions in the scalar potential during inflation. The subsequent 
coherent oscillation of the field obtains a phase rotational motion and begins to carry 
non-zero baryon number in the presence of 5-violating operators. A crucial point is 
that this coherent oscillation of the field is unstable, with spatial perturbations and 
fragments into the Q-balls after dozens of oscillations. Recent detailed lattice simulations 
have revealed that almost all the baryon asymmetry initially carried by the coherent 
field is absorbed into Q-balls || 0- 

The large amplitude of the field inside the Q-ball protects it from being thermalized. 
The typical decay temperature of Q-balls is given by T d ^(a few) GeV, which is well below 
the freeze-out temperature of the LSP. As a result, the LSPs produced by the Q-ball decay 
are not thermalized and retain the initial abundance. Therefore, the mass density of the 
non-thermal neutralino dark matter Q x is directly connected to that of baryons : 



where m x and m p are the LSP and nucleon masses, respectively; ng and are the baryon 
and field number densities. One can easily understand this relation by considering that 
the decay of a single field produces ^ 1 neutralino LSP in the final state. Here, the ratio 
of ub to is fixed when the field starts coherently oscillating and remains constant 
until it eventually decays. Under the assumption of i?-parity conservation, the possible 
maximum value of this ratio is (ns/n^) = 1/3. The typical value of this ratio, which 
naturally appears in models of AD baryogenesis in the mSUGRA scenario is ns/n^ ^ 0.1. 

Unfortunately, this simple model is not cosmologically viable. We need an extremely 
light bino m x ^ 1 GeV to explain the required mass density of dark matter. This fact 
indicates that the formation of a Q-ball is a serious obstacle for the AD baryogenesis 
rather than a solution to the Q B -Q DM coincidence puzzle. 

There are several ways to reconcile this situation. The obvious one is to make the 

Q-balls small enough that they can evaporate well above the freeze-out temperature of 

the LSP. This can be done by gauging the XJ(l) b-l symmetry [|7| or assuming a relatively 

strong three-point coupling of the inflaton to the gluino ||. In both cases, it is clear that 

the relation between fig and fl x is lost completely. Another possibility is to adopt the LSP 

2 In the original work, the baryon absorption into the produced Q-balls was assumed not to be so 
efficient, which conflicts with the detailed lattice simulations @, n|. 
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with a large annihilation cross section, such as higgsino H or wino W. In Refs. ||, pTCf] , 
it is shown that the subsequent pair annihilations of these LSPs can naturally lead to 
the desired mass density of dark matter. [] However, in this case, one expects that the 
resultant mass density of baryons and that of dark matter are completely independent 
of each other, since the number of produced LSPs should be drastically reduced via 
annihilation processes so as not to overdose the Universe. However, astonishingly, this is 
not always the case. 

In this letter, we will show that a class of AD baryogenesis scenarios allow us to 
connect the two crucial quantities in our Universe. Our final goal is to show the following 
relation between the mass density of baryons and that of dark matter: [] 

n x VMx / \ m x) 

where is the mass of the flat direction field <p, which is approximately given by the 
squark mass; (av ) x is the s-wave component of the annihilation cross section of the LSP; 
5 c g = 0(0.1) is the effective CP-violating phase of the <p field. A numerical factor in 
the right-hand side depends on the details of the produced Q-balls, but basically they 
are calculable. If we take the typical annihilation cross section of the H- or HMike LSP, 
(crv) x ~ 10~( 7 ~ 8 ' GeV~ 2 , one can easily see that this is just the desired relation. In the 
remainder of this paper, after a brief review of the AD baryogenesis, we will derive this 
interesting relation and discuss the conditions for it to hold. 



2 The model 

Our model of the AD baryogenesis is basically the same as the one presented in Section 
II B of Ref. For self-consistency, we briefly review the mechanism with particular 



attention to the late-time decays of Q-balls. The derivation of the relation Eq. (|2]) and 
some conditions for it to hold will be discussed in subsequent sections. 

We consider the situation that there exist some chiral symmetries, such as i?-symmetry, 

3 The LSP is not necessarily a pure H or a pure W . A significant mixing of the bino component in 
the LSP is quite possible. In the case of a large tan/3, the acceptable mass density of dark matter can be 
obtained even with the LSP whose dominant component is a bino. 

4 As we will see, when the mass of the gravitino to 3 / 2 is larger than m^, in Eq. (^) should be 
replaced by m.3/2. 
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which forbid non-renormalizable operators in the superpotential that lift the relevant flat 
directions for baryogenesis. B- violating operators needed for AD baryogenesis are supplied 
in the Kahler potential. We take the following operators, for example, which are consistent 
with the i?-symmetry: 



(3) 



where M* = 2.4 x 10 18 GeV is the reduced Planck scale, and Ai,A2 are the coupling 
constants of the order of 1, which are generally complex numbers. Here, Q, U, D, L, I 
and Z denote superfields (and their scalar components) of left-handed quark doublets, 
right-handed up-type and down-type quarks, left-handed lepton doublets, inflaton and 
the one relevant to SUSY breaking in the true vacuum, respectively. These interactions 
induce the following terms in the scalar potential: 

/ ott2 3 m 2 \ 

SV=- U—QW&L + \ 2 ^QUW^L + h.c.J , (4) 

where is the mass of the gravitino, and H is the Hubble parameter of the expanding 
Universe. 

If the expectation value (QWD^L) is large enough after the inflation, the terms in 
Eq. (|4]) give the phase rotational motion to this scalar condensate and generate net baryon 
asymmetry. In fact, in the scalar potential of the minimal supersymmetric standard model 
(MSSM), there are many D- and F-flat directions that give a non-zero expectation 



value of (QU' D< L) . In the remaining discussion, we adopt the flat direction labelled by a 
linear combination of the monomials of the chiral superfields: UDD, QDL. An explicit 
parametrization of this direction, as suggested in the original work H, is given by 




= (<!>* Ok EjM = ( <fc 0], = ( 2 ) , 

(5) 

where superscripts (l)-(3) denote the generations of the chiral superfields. Here, the 
column and low vectors denote the isospin of the SU(2) L and the colour of the SU(3)c 
gauge groups of the MSSM, respectively. Hereafter, we use this simple parametrization. 
5 The £>-violating operator oc QQU^E^ can also be applied to the following arguments. 
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Although there are many other parametrizations of the flat direction, the main arguments 
in the following discussion do not change much. 

To give large expectation values to the flat-direction fields, we assume the four-point 
couplings of these fields to the inflaton in the Kahler potential |12|] : 

6K =TH + & 20202 + &3030 3 ) , (6) 

where fej's are real coupling constants of the order of 1. Here we use the same symbols for 
the chiral superfields as for the corresponding scalar components in Eq. (|5|). By calculating 
the scalar potential of the supergravity, we can find that the following terms are induced 
at leading order: 

5V = 3[1- (h + b 3 )} tf 2 |0i| 2 + 3 [1 - (h + 6 3 )] # 2 |0 2 | 2 . (7) 

Therefore, if both {b\ + 63) and (62 + ^3) are somewhat larger than 1, 0i and 02 (and so 
(f) 3 ) obtain large expectation values during inflation. Furthermore, if (b\ +b 3 ) « (6 2 +63), 
we can naturally expect that \<pi\ ~ \<f>2\ (= |0|). We assume that this is the case in the 
remainder of this paper. 

Now, we can approximately write down the relevant scalar potential of the </> field as 
follows: 

V = K - chH 2 )^ + ~^ 2 M 4 + h.c.) + ^||K0 4 + h.c.) + . . . , (8) 

where the ellipsis denotes the higher-order terms coming from the Kahler potential; a#, 
a m are complex coupling constants; c# = 0(1) is a real coefficient; denotes the soft 
SUSY-breaking mass of the <f) field, which is roughly given by the average value of the soft 
SUSY-breaking masses of the Q, U, D and L. Here, we have neglected the terms induced 
by the interactions with thermal backgrounds, which will be justified later. 

We are, now, at the point where we want to discuss the evolution of the <fi field. During 
inflation, the large negative Hubble-induced mass term causes an instability around the 
origin, and the (ft field develops a large expectation value. The amplitude of the field just 
after the inflation, |0 O |, is determined by the balance point between the Hubble-induced 
mass term and non-renormalizable operators in the Kahler potential. In the following 
discussion, we treat it as a free parameter and assume l^lo^M*. (Actually, there exists 
an interesting method to fix |</>|o below the Planck scale. This can be done by gauging 



5 



the U(1)s_l symmetry. In this case, we can fix the amplitude of the field at the B — L 
breaking scale, |0|o — vb-l- Q The details on this point are discussed in Ref. 0.) 

After the end of inflation, the amplitude of the Hubble parameter gradually decreases. 
When exceeds the Hubble parameter, the field starts coherent oscillation around 
the origin. At this time H = H osc ~ m^, a huge baryon asymmetry is produced because 
of the second and third operators in Eq. (||). As long as |0|oi$M*, the curvature along 
the phase direction is smaller than the Hubble parameter during the inflation. Therefore, 
the initial phase of the field is generally displaced from the bottom of the valley of the 
scalar potential. This is the reason why the second and third operators give the phase 
rotational motion to the field. 

It is not difficult to estimate the ratio of baryon to 0-number density, which is fixed 
at H = Hqsc , as 




max 



m[ m:) W UJ U* 1 eff 



(9) 



where 5^ s = su^arga^ + 4arg(0 o )), and 5^ = sin(arga m + 4arg0 o ). [] Note that this ratio 
remains constant until the field eventually decays into the SM particles and the LSPs. 

After dozens of oscillations, the scalar condensate of the field fragments into the 
Q-balls. Almost all the baryonic charge and energy density carried by the coherent field 
are absorbed into the Q-balls [|, As for the details about the Q-ball, see Ref. [|ItJ. 
The produced Q-balls behave as ordinary matter, and their energy density decreases as 
Pq(= p^) oc R~ 3 , where R is the scale factor of the expanding Universe. On the other 
hand, the energy of the inflaton is converted into radiation through reheating. After 
completion of the reheating process, the energy density of the radiation decreases as 
Pr oc R~ a . Then, the energy density of the Universe is dominated by the produced Q- 
balls before their decays, which take place at T = T^, when the following condition is 
satisfied: 

T « >3T *(wS' (10) 



6 Although the operator conserves the B — L symmetry, the flat directions carry non-zero B — L charges 
with the same sign; then, they can be lifted at the B — L breaking scale by the potential induced by the 
U(1) B -L £>-term §. 

7 This ratio cannot exceed 1/6 for our choice of the flat direction. If m 3 / 2 3> m^, as in the case of 
anomaly-mediated SUSY-breaking models, the condition |</>|o ^ m^M^/m^^ is necessary for the <f> field 
not to be trapped at the local (or global) minimum located near the Planck scale Jll| . 
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where Tr is the reheating temperature of the inflation. Although it depends on the 
size of the Q-ball charge, the typical decay temperature of the Q-ball lies in the range 



10 MeV^Td^(a few) GeV, and hence the relation Eq. ( |T0|) can be easily satisfied. In 




the following, we assume this to be the case. We will come back to this condition after 
the derivation of Eq. (§). 

The resultant baryon asymmetry after the decays of the Q-balls is now given by the 
following simple formula: 

n B 
s 

where s is the entropy density of the Universe. Here, in the second equality, we have 
used the fact that the effective mass of the Q-ball per 0-number is given by in a good 
approximation. Note that the resultant baryon asymmetry is completely independent of 
the detailed history of the Universe, since the Q-balls dominate the energy density of the 
Universe before they decay. 

If the LSP x nas a large s-wave annihilation cross section, which is the case for the 
H- or VT-like LSP, f\ the resultant relic abundance of x is given by the following simple 
expression to quite good accuracy []9], |10|] : 



45 



{av) 



-i 
x 



8n 2 g,(T d ) M*T d 



(11) 



where n x is the number density of the LSP, and g*(Td) denotes the relativistic degrees of 
freedom at T = T d . In terms of the density parameter, the resultant matter density of 
baryons and that of dark matter can be expressed as follows: 



Q B h 2 



0.02 




V 




rG 



x 10 fc 



(12) 




^100 MeV^j AO" 7 GeV- 2 * 



(av) 



(13) 



x 



8 The sneutrino v is also an interesting candidate. Even the bino-like LSP can have a large s-wave 
component in the case of very large tan/? via H contamination. 
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Therefore, if either H or W is the LSP, the baryon asymmetry and the dark matter may 
be explained simultaneously in the following range of parameters: 

10 15 GeV < |4, < 10 17 GeV , 10 MeV < T d <(a few) GeV , 

10- 8 GeV- 2 <(av) x < 10- 7 GeV" 2 , 10" 2 < \a\6 eS < 10" 1 , (14) 

where a and S e s denote either a# or a m , and either 5^ or S^g, respectively. 



3 Derivation of the relation of Qb an d Qdm 

Now, we come to deriving the relation in Eq. (|2]). For this purpose, we have to know the 
|0|o dependence on the decay temperature T d of the Q-ball. The terms \(p\ and T d are 
related to each other through the size of the Q-ball charge U Q" . 

The size of the Q-ball charge crucially depends on the initial amplitude of the field 
and its scalar potential. The scalar potential relevant at the time of the Q-ball formation 
can be written as 

W) = ^(l + ^log(^))|0| 2 , (15) 

where Mq is the renormalization scale at which the soft mass is defined, and the 
inog(|0| 2 ) term represents the one-loop correction. This mainly comes from the gluino 
loops and the typical value of K is found to be in the range -0.1 <K< -0.01 JLJ. This 
negativeness of the i^-factor makes the potential a little bit flatter than the quadratic one, 
which is a necessary and sufficient condition for the Q-ball to be formed in the present 
baryogenesis. 

Recently, the typical size of the Q-ball charge has been calculated by detailed lattice 
simulations [|], ^j. Applying their result to our model, the size of the Q-ball charge can 
be estimated to 

Q=^y e , do) 

\ m 4>) 

where j3 ~ 6 x 10 -3 , which depends on the .fT-factor and the fluctuations of the field 
after inflation, and e ~ 0.01 is a constant independent of \4>\q. F\ 



9 This is the case for (ns/n^) 0.01, which is naturally satisfied in the present model with |</>|o M* 
If this is not the case, e should be replaced by (ns/n^). 
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The size of the charge that evaporates from the surface of a single Q-ball through 
interactions with thermal backgrounds is about AQ ~ 10 18 [|1(J. Therefore, as long as 
|0|o~ 10 14 GeV, the Q-ball survives thermal evaporation. The remaining charges of the 
Q-ball are emitted through its decay into light fermions. The decay rate was calculated 

u 3 A 



in Ref. |17| as 



To = ?- 



dQ < 

dt ~ 192tt 2 



(17) 



Here, A = 4:TiRq is the surface area of the Q-ball, where Rq ~ \/2/ {m^J\K\) is the Q- 
ball radius; u ~ is the effective mass of the Q-ball per 0-number. This upper bound 
is likely to be saturated for 0(0) ^> m^, where 0(0) is the field value of at the centre of 
the Q-ball. This is the case in the present model. 

Then, the decay temperature of the Q-ball is given by 




1/4 fm,MA 1/2 



1 TeV 



V2 /10 



20' 



1/2 



Q 



(18) 
(19) 



where rj ^ 1 denotes the ambiguity coming from an inequality of the decay rate in Eq. (|T^) . 
Then, in terms of the initial amplitude of the field, the decay temperature is written as 
follows (see Eq. flU)): 



90 



, 1/4 /m,M*\ 1/2 



487^1/36 \^ 2 9*{ T d 



(20) 



From Eqs. (|T2"D, ( pT3[ ) and (p0|), we can see that both the matter density of baryons 
and that of dark matter are linearly proportional to |</>|o, and hence we can easily derive 
the wanted relation: 




16n\K\Pe 



x max < 





a H \S§{ 




< 

\a m \ 


( m 3/2 \ 
V m 4> J 


2 



(21) 



10 This factor is, in principle, numerically calculable by determining the accurate profile of the Q-ball, 
which requires the details of the scalar potential. There might also exist an ambiguity from the decay 
channels induced by loop diagrams. 



9 



If we assume (\an\, \a m \, v) ~ 1 an d the typical values for the parameters, K, (3 and e, 
the above equation is written as 



^£ ~ 10 3 - 4 ( ay 

\ m x / 



x max 



m l \ sH ( m 3/2 \ s g 
°eff J /_„.\-1 °eS 



M-v ett ' H; 1 



(22) 



By using the typical annihilation cross section for H- or VT-like LSP, (crv) x ~ 5 x 
10~ 8 GeV~ 2 , we derive 




5 x 10" 8 GeV~ 2 




(23) 



where m and 5 e s denote either or m 3 / 2 , and either 6*g or S^ s , respectively This is 
perfectly consistent with observations. 

We can obtain another interesting piece of information by "multiplying" the mass 
density of baryons and that of dark matter together. In this case, we can remove the 
ambiguity associated with the decay temperature of the Q-ball. It is easy to show that 
the following relation holds: 



n B h 2 x n x h 2 = 1.76 x 10- 2 ^ . (24) 




This relation gives us information on |0| o , the initial amplitude of the <fi field. If we take 
the typical values for parameters, we can see that the present mass density of baryons 
and that of dark matter suggest 10 15 GeV ^ \4>\q ^ 10 16 GeV, which surprisingly coincides 
with the B — L breaking scale suggested from the see-saw neutrino masses |Tj|. Note that 
this value does not affect the prediction in Eq. fl2"ip. 



4 Conditions for the relation to hold 

In this section, we discuss the conditions for relation (|21~D to hold. We have to examine 
the thermal effects, on the scalar potential, of the field given in Eq. (|J). If the thermal 
effects dominate the scalar potential when the <fi field starts to oscillate, they cause the 
evaporation of the (f) field before the formation of the Q-balls, or drastically change the 
|</>|o dependence on the Q-ball charge. 

First, let us discuss the effects related with the thermal mass terms. If the cosmic 
temperature T satisfies f\<p\o < T, the field coupled to the cf> field through the coupling 
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constant / induces the thermal mass term c 2 / 2 T 2 |0| 2 , where c is a real constant of the 
order of 1. Therefore, if the two conditions, f\<fi\o < T and cfT > H, are satisfied 
simultaneously when H ^ m^, the thermal mass term causes the early oscillation of the </> 
field |L2|) |ig| , [20[| . In this case, the field is thermalized and Q-balls are not formed. The 
condition for the reheating temperature of inflation, Tr, to avoid this early oscillation is 
given by 

This condition should be satisfied in all the gauge and Yukawa coupling constants asso- 
ciated with the field. The most stringent constraint comes from the first term with 
Yukawa coupling of an up quark, Y u ~ 10 -5 . This constraint becomes significantly weak 
if we adopt the second and third generations of squarks for the flat direction. 

A much more stringent constraint comes from the thermal logarithmic potential |H 




5V D a g a T 4 log ™ , (26) 



where \a g \ = 0(1), and a is a constant given by the fourth power of gauge and/or Yukawa 
coupling constants. This leads to the following constraint on the reheating temperature J7| 
0: 

r R <^^W (27) 



a g \a 



From Eqs. ( pq) and ([27]) combined with the Q-ball dominance condition in Eq. (|T0|), 
we obtain the allowed region of the reheating temperature where the relation Eq. (f2~l|) 
holds. We show this region in Fig. [I]. 

From this figure, we see that there is a wide allowed region for the reheating tem- 
perature if |0|o~ 10 15 GeV. Note that, even if we use Tr 2> 10 8 GeV, there is no "cos- 
mological gravitino problem" [23[] as emphasized in Ref. [10|. This is because the large 



entropy production associated with the decays of Q-balls dilutes the gravitino number 
density substantially. As long as the reheating temperature and the initial amplitude of 
the field are within this allowed region, Qb/^x * s fi xe d by the low-energy parameters, 
and is totally independent of Tr and \4>\q- 



11 Here, we have used the fact that the cosmic temperature behaves as T ~ (iJT^Af*) 1 / 4 before the 
completion of the reheating process of inflation. 
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Figure 1: The allowed region of reheating temperature for relation ([H]) to hold. The 

= 1/20. 



red (solid) line denotes the upper bound from Eq. (|2~7|) with \a g \ = 1 and 
The purple (dotted) line denotes the upper bound from the thermal mass term, Eq. (|25|) . 
Here, we take c = 1, / = 10~ 5 . This constraint becomes significantly weak if we do not 
adopt the first generation of up-type squark in the flat direction. The blue (dashed) lines 
represent the lower bounds from the Q-ball dominance condition given in Eq. (ffD) with 
r\ = 0.3, \K\ = 0.1 and 77 = 1 , \K\ = 0.03 from the bottom up, respectively. Throughout 
this calculation, we have used = 1 TeV. 

5 Conclusions and discussion 



In this paper, we have proposed a solution to the VL B -VL DM coincidence puzzle. We have 
pointed out that a class of the AD baryogenesis directly relates the observed mass density 
of baryons to that of dark matter with the help of the late-time decays of Q-balls, if the 
LSP has a large s-wave annihilation cross section, as in the case of H- or IV-like LSPs. 
The relation we have shown is 

(28) 

A beautiful point of this scenario is that the ratio is totally independent of the reheating 
temperature of inflation and the initial amplitude of the <fi field, as well as from a detailed 
12 In the case of m 3 / 2 > m^, we should replace by m 3 / 2 - 
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history of the Universe, once the relatively weak constraints shown in Fig. |I| are satisfied. 
We should stress that the ratio is solely determined by the low-energy parameters, except 
for the O(0. 1) effective CP phase, which makes this scenario quite testable in the future 
experiments. Our scenario is free from the "cosmological gravitino problem", since the 
large entropy production associated with the Q-ball decay sufficiently dilutes the number 
density of gravitinos. 

In addition to the above beautiful relation, there are many interesting implications of 
the present scenario. The large H or W component of the LSP significantly enhances its 
detection rates in both direct and indirect dark matter searches [p!0[ . Furthermore, such 
non-thermal dark matter may naturally explain the observed excess of positron flux in 
cosmic rays f24|. 
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